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Abstract 

It is a matter of experience that nonlinear waves in dispersive media, propagating primarily 
in one direction, may appear periodic in small space and time scales, but their characteristics — 
amplitude, phase, wave number, etc. — slowly vary in large space and time scales. In the I970’s, 
Whitham developed an asymptotic (WKB) method to study the effects of small “modulations” 
on nonlinear periodic wave trains. Since then, there has been a great deal of work aiming at 
rigorously justifying the predictions from Whitham’s formal theory. We discuss recent advances 
in the mathematical understanding of the dynamics, in particular, the instability of slowly 
modulated wave trains for nonlinear dispersive equations of KdV type. 


1 Introduction 

The modulational instability, at the core, is the observation that nonlinear wave trains in a dis¬ 
persive medium may be unstable to self modulation, developing nontrivial large-scale structures, 
which continue to evolve as they propagate. The phenomenon is familiar to researchers in many 
fields. Moreover, it is a fundamental issue in the theory of nonlinear waves. The history of the 
modulational instability of nonlinear dispersive waves is long and fascinating. We do not have 
space here for a complete account of it, and we will discuss only selected high points. For more 
details, the interested reader is encouraged to read the excellent review article by Zakharov and 
Ostrovsky m- 

The first understanding of the importance of the modulational instability arouse in hydrody¬ 
namics in the 1960’s. Whitham m applied his newly developed theory of nonlinear dispersion to 
Stokes’ periodic waves in water, associating the hyperbolic (or elliptic) nature of the resulting mod¬ 
ulation equations with the stability (or instability, respectively) of the underlying wave; he showed 
that a Stokes wave is unstable if kho > 1.363 ... , where k is the wave number of the underlying wave 
and /iQ is the undisturbed fluid depth. At about the same time, Benjamin and Feir [7], motivated 
by laboratory experiments [26], derived their eponymous instability; they showed that a Stokes 
wave in deep water is unstable. Additionally they discovered the sideband nature of the instability 
— the marginal stability of the wave number, with a band of unstable wave numbers to either 
side. Corroborating results arrived nearly simultaneously, but independently, by Lighthill [57], 
Ostrovsky |60j . Benney and Newell [8], and Zakharov [zaizi]. 

The modulational instability has been experimentally observed, either directly or indirectly, in 
many physical systems. For example, pulses in nonlinear optical fibers with so-called anomalous 
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dispersion are (approximately) governed by the focusing, cubic nonlinear Schrodinger equation, 
which exhibits the modulational instability. The importance of the effect was noted in early papers 
by Anderson and Lisak [2] and Hasegawa m. where it was proposed as a method for generating 
short optical pulses; it was later observed experimentally by Tai, Hasegawa and Tomita [67]. Since 
then, there have been a great deal of work observing the phenomenon in a number of optical 
settings; see [55], [68], [30], [15], [65], for instance. 

The modulational instability and associated Whitham equations for integrable systems, such as 
the KdV and nonlinear Schrodinger equations, were studied by Flashka, Forest and McLaughlin |28j 
and by Lax and Levermore [56], and it has grown into a large and active field of research; see [25], 
[25], [22] . [54], [69], [33], for instance. Moreover, there has been a great deal of recent work aiming at 
making the connection between the elliptic nature of the modulational equations and the instability 
of the underlying wave. To single out a few in the vein, we mention m, m, [m, m- Here we 
will be mainly concerned with the latter program, how Whitham’s formal theory translates into 
rigorous stability results. This has led to many applications in nonlinear dispersive equations; see 
[35], m, [3l], m, m, Ha, m, M, m, Hl], for instance. 

We begin in Section 2 by constructing periodic traveling waves to the generalized KdV equa¬ 
tion. In Section 3 we will review Whitham’s theory of modulations and discuss implications for the 
modulational instability. Section 4 contains the main analytical contributions, providing a rigorous 
mathematical program to determine the modulational instability of periodic traveling waves of KdV 
type equations. For illustrative purposes, in Section 5, we will discuss applications to the KdV, 
modified KdV and Schamel equations. We will then extend the theory to equations admitting non¬ 
local dispersion and discuss applications to the Benjamin-Ono equation, as well as small amplitude 
waves for the fractional KdV, intermediate long-wave equations and Whitham’s equation for water 
waves. 

2 Periodic traveling waves of generalized KdV equations 

Much of the chapter is devoted to the modulational instability of periodic traveling waves to the 
generalized Korteweg-de Vries (KdV) equation 

(2.1) Ut — Uxxx T /{ujx- 

Here t G M is typically proportional to elapsed time and x G M is related to the spatial variable in 
the predominant direction of wave propagation; u = u{x, t) is real valued, representing the wave 
profile or a velocity, and / is a suitable nonlinearity. In many examples of interest, / obeys a power 
law. Throughout we express partial differentiation either by a subscript or using the symbol d. 

Perhaps the best known among equations of the form (12.111 is the KdV equation 

Ut — Uxxx T {u 

itself which was put forward in m and [53] to model the unidirectional propagation of surface water 
waves with small amplitudes and long wavelengths in a channel. It has since found relevances in 
other situations such as Fermi-Pasta-Ulam lattices; see m, for instance. In the case of f{u) = 
au^-\-bu^, a, b constants, (|2.ip is called the modified KdV or Gardner equation, which models internal 
waves propagating in a density stratified fluid; see [M], for instance. In the case of f{u) = 
moreover, (I2.1|l is called the Schamel equation; see [M], for instance. 

By a traveling wave, we mean a solution which propagates at a constant speed without changing 
the shape. Therefore it takes the form u{x,t) = u{z), z = x — ct, where c G M is the wave speed 
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and n is a t-independent solution of 


Ut = U^zz + CU^ + f{u)z- 


It then reduces by quadrature to 

(2.2) + CU + f{u) = a 


for some constant a G M. It further reduces by quadrature to 


(2.3) 



E + au - cu^ 

2 


F(u), 


where F' = f, F{0) = 0, for some constant i? € M. We refer to ()2.3I) as the profile ODE cor¬ 
responding to (j 2 .I|) . One may employ elementary phase plane analysis to infer the existence of 
periodic, homoclinic (pulses) and heteroclinic (fronts) orbits of (|2.3p . Note that the existence and 
non-existence of bounded solutions of ( 12 .3|) is determined by the effective potential 


V{u; a, c) 


+ -cu^ 


— au. 


2.1 Some Explicit Solutions 


In some special cases, (12.3p may be directly integrated to yield classes of explicit solutions. In the 
case of / = 0 in ( 12 . 11 ) . for example, namely the (linear) Airy equation, ( 12 .ip admits plane wave 
solutions u{x,t) = where k represents the wave number and c is the wave speed, provided 

that the temporal frequency w := fee is related to fe via the dispersion relation uj = — fe^. We will 
return to such linear waves in Section 3.1 when we discuss Whitham’s theory of modulations. The 
modulational stability of these waves, it turns out, is trivial and follows directly from the dispersion 
relation. In fact, all these waves are modulationally stable. Therefore we will be concerned with a 
genuine nonlinearity, typically f{u) = for some p > 0 and a ^ 0 real. Thanks to the scaling 

invariance of m, by the way, a may be taken arbitrary. In fact, the stability of periodic traveling 
waves of (j 2 .ip is independent of the choice of a. 

In the case of f{u) = namely the KdV equation, moreover, simple phase plane analysis 
implies that periodic solutions of (|2.3p exist, depending on a, E, and c, provided that the cubic 
polynomial E — V{■] a, c) has three distinct real roots, which we label as 7 < /3 < a. Periodic 
solutions then correspond to oscillation between in the interval [a,/3], and o, E, and c are related 
to a, /3, and 7 via 


(2.4) 


E = 


al3'y 


a = -- 

,2 


aP + (3'^ + aj 

~~6 ’ 


c = 


a + (3 + 'y 


In particular, ()2.3I) may be written as = {a — u){u — I3){u — 'j)/3 and u may be defined implicitly 
via 

z + zo dw 


Vs Ju -^(a — w){w — f3){w — 7 ) ’ 

where zq is a constant of integration associated with translational invariance and can be chosen so 
that tt(0) = a. We then make the change of variables 


w = a — {a — (3) sin^ 


sm0 = 


a — w 
w-V 
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and solve the above for u explicitly as 


(2.5) u{z;a,/3,'y,zo) = P + {a-/3)cn^(^ 

where cn(-; m) denotes the Jacobi cnoidal function with elliptic modulus m E [0,1). These cnoidal 
waves define a four-parameter family of periodic solutions of the KdV equation. Note that the 
period is T = (m) /\/a — (5, where K is the complete elliptic integral of the first kindil- 

Another example of explicit periodic traveling waves is when f{u) = in (|2.1I) . which are 

called the focusing (+) and defocusing (-) modified KdV equations. In the focusing case, V{- ;a, c) 
describes a bistable double-well potential and hence simple phase plane analysis implies that non¬ 
constant periodic solutions of (j2.3l) exist, depending on a, E, and c, provided that the quartic 
polynomial E — V{-; a, c) has four distinct roots (i.e. its discriminant is non-zero), say, a, 13, 7 
and 6, with at least two of them being real; see Figure [H for instance. The general solution is 
rather tedious with a number of cases to consider, but it amounts to reducing the elliptic integral 
to standard form. In fact, if we write the solution in the form 

du 

— ,, - = uit 

y/{u- a){u - j3){u - 7 )(m - S) 



for some appropriate constant a; E M and make the Mobius tranformation 


au + b 
V = - 

u — a 


for a, b constants, which maps one of the roots of the quartic to 00 , then the above elliptic integral 
becomes 

f dv 

/ — / = 

J 

for some cubic polynomial P 3 , which may be analyzed as above. Clearly, the root a must not be 
an end point of the interval over which the solution oscillates. 

In the case of a = 0, in particular, the effective potential V{- ;0, c) is symmetric (even) and the 
expressions for explicit solutions simplify greatly. Specifically, there are two families of solutions, 
one in terms of the Jacobi elliptic cnoidal function cn and the other in terms of the Jacobi dnoidal 
function dn. The cnoidal waves correspond to solutions with E > 0 and their orbits in phase space 
lie outside the separatrix; see Figure [T](a). The dnoidal waves correspond to solutions with E < 0 
and their orbits in phase space are bounded by a homoclinic orbit (corresponding to a solitary 
wave). In the case of c < 0 and E < 0, for instance, the quartic on the right side of the profile ODE 


9 I—? 9 ^ 

= 2E — cu^ - 

^ 6 


has four real roots, u = ±fci,u = 3zk2, with ki = —3(c + y + ^) and k 2 = —3(c — y + ^)- 
The change of variables u = k 2 V gives 

^ Note that cn( 2 :;m) —>■ sechfz) and K{m) —^ 00 as m —^ 1“ and hence (12.51) reduces to the well-known formula 
for the KdV soliton. 
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Figure 1: Representative graphics of the effective potential energy V{- ]a,c) in the case of the (a) 
focusing and (b) defocusing modified KdV equations. In the focusing case, periodic solutions exist 
only for E > E 2 {a, c). There are two families of periodic solutions when E G (£' 2 , F^i), while there is 
only one family when E > Ei. In the defocusing case, periodic solutions exist only for E G (E 2 , £ 1 ) 
and there is only one family of solutions. 


Note that the elliptic function y{z) = dn( 2 ;; k) solves 

Therefore v = dn{^z] k), = 1 — The other case, when c < 0 and £ > 0, which is solvable 

in terms of the cnoidal function, follows similarly. 

In the defocusing case, non-constant periodic solutions exist, depending on a, £, and c, only 
when the quartic polynomial E — V{- ;a,c) has four distinct real roots; see Figure [Hb). In the case 
a = 0, the potential is likewise symmetric and one may find explicit solutions in terms of the Jacobi 
snoidal function. 

2.2 General Existence Theory 

While an extensive amount of literature is devoted to studying properties of explicit, periodic trav¬ 
eling waves of m, e.g. the elliptic function solutions of the KdV and modified KdV equations 
discussed above, our modulational instability theory does not make use of such formulae. (Never¬ 
theless, we will strive to state our general results in the context of explicit solutions when possible.) 
Rather, of importance to us are properties of the parameterization of periodic traveling waves with 
respect to various coordinate systems and the Jacobians between them. As a first step in the direc¬ 
tion, we remark that the Jacobian of the map (12.4p from to itself is non-singular, provided that 
the cubic polynomial E — V{-;a,c) has three distinct roots, a property fundamentally related to the 
sign of its discriminant. As a result, in the set where the discriminant of the polynomial is positive, 
£ — V(-;a, c) has three distinct real roots and the transformation (12.4p be smoothly invertible. 
Consequently, the cnoidal wave solutions (12.5p of the KdV can be considered as a four-parameter 
family of periodic traveling waves parametrized by the “natural” ODE parameters a, £, c, and zq, 
whose profile and fundamental period depends smoothly on these parameters. Analogous results 
hold for elliptic function solutions of the focusing and defocusing modified KdV equations: within 
the domain of existence they can be smoothly parametrized in terms of £, c, and zq (recalling that 
a = 0 in these explicit solutions). 

With the above in mind, we note that if / in (12.11) is a polynomial nonlinearity then the zero 
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Figure 2 : Schematic drawing showing the locations of the roots 


and . 


set of the discriminant oi E — !/(•; a, c), denoted as 

r := I (a, Fi, c) G : disc + au — — F{u^ ~ ’ 

defines a variety dividing the parameter space into open sets with a constant number of periodic 
solutions. Solutions of ()2.3p with (a, E,c) € F correspond to either equilibrium or homoconic 
solutions of ()2.3p . which, respectively, correspond to constant solutions and solitary waves of (12.11) . 
In particular, (0,0, c) G F represents the solitary waves homoclinic to zero. 

In order to ensure the existence of periodic orbits of (12.3p . we take {a,E,c) G \ F so that 
there exist simple and distinct roots of the equation E = V{-\a,c) and that there are real roots u± 
satisfying u_ < u+ such that V{u‘,a,c) < E foi u G (ri_,tt+); see Figure[2l It follows that the roots 
u± are smooth functions of the traveling wave parameters (a, E, c) and, without loss of generality, 
we can assume n(0) = U-. It follows that the corresponding periodic solution of (12.31) is an even 
function of z = x — ct and that its period can be expressed via quadrature as 


T = r(a, E, c 



du 

y/E -V{u-,a,c) 


y/2 / du 


where integration over 7 represents a complete integration from u- to and then back to U- 
again. Notice however that the branch of the square root must be chosen appropriately in each 
direction. Alternatively, the contour 7 can be interpreted as a closed loop (Jordan curve) in the 
complex plane that encloses a bounded set containing both u±. By a standard procedure, the above 
integral can be regularized at the square root branch points and hence represents a function of 
a, E, and c. 

In general, (12.ip admits three conserved quantities, physically corresponding to the Hamiltonian, 
the mass, and the momentum. Given a T-periodic traveling wave u{--,a,E,c) of (El]), they are 
defined by 


( 2 . 6 ) 


H{a, E, c) 
M{a, E, c) 
P(a,E,c) 


(!i _ f („)) = k! / r-V(n-,a,c)-F(u) 

V 2 ^ V 2 L JE-V(irar\ 


f 

f 


u dx = 


VI I _ 

2 ^yE — V{u; a, c)' 

2 , \/2 / du 

u dx = — ® , 

2 V y^E-V{u;a,c) 


V y/E - V{u] a, c) 

u du 
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respectively, where the integral over 7 is defined as before. As above, these integrals can be 
regularized at their square root branch points and hence represent functions of the traveling 
wave parameters. As we will see, the gradients of the period and the above conserved quantities 
plan a fundamental role in the stability of the periodic traveling wave solutions of () 2 . 1 [) . 

In summary, (12.111 in the case of polynomial nonlinearity admits a four-parameter family of 
periodic traveling waves, smoothly parametrized by a, E, c in (I2.3jl and an auxiliary parameter zq 
associated with the translational invariance of the equation. Solitary waves, whose profile asymp¬ 
totically vanishes, correspond to a = E = 0, and hence they form a co-dimension two subset of the 
full set of traveling waves. 

3 Formal Asymtotics and Whitham’s Theory of Modulations 

We are interested in the stability of periodic traveling waves of nonlinear dispersive equations 
to slow modulations. That means, we wish to understand the long-term dynamics of “slightly” 
modulated periodic traveling waves, in particular, to determine whether the long-term dynamics 
will be appropriately described by small modulations of the carrier wave. In this section, we 
illustrate some asymptotic calculations that shed light on the long-term dynamics and hence on the 
modulational stability. Mainly, we will utilize a formal multiple-scales expansion to derive a first 
order system of equations governing the slow dynamics of averaged quantities of the modulated wave 
over the fast scale. Such an approach was pioneered by Whitham in the 1960’s and early 1970’s 
(see [70], [TT], [72], for instance), and it has since been widely applied in both the physics and 
mathematics communities. As such, the formal theory has been affectionately termed Whitham’s 
theory of modulations, or simply, Whitham’s theory. 

To illustrate the key ideas, we begin by considering slow modulations of a linear dispersive 
wave packet, where the modulational dynamics are simplified due to the direct coupling of spatial 
and temporal frequencies via the linear dispersion relation. We will then extend this calculation 
to nonlinear wave packets where one must make use of a nonlinear dispersion relation to close the 
modulation system. Section 14.3[ we will discuss rigorously validating the predictions from these 
calculations. 

3.1 Linear Dispersive Waves 

Notice that any constant function u{x,t) = uq is a solution of (I2.I|1 . Seeking nearby solutions of 
the form 

u{x,t) = uq + ev{x,t), |e| <C 1, v = 0(l), 

we arrive at that 

T '^xxx T f (^o)w — 

Ignoring all 0(e), we then find that as long as v remains 0(1) we can approximate u by a solution 
of the linear equation 

(3.1) Vt + Vxxx + f(uo)vx = 0. 

It admits plane wave solutions of the form 

(3.2) v{x,t) = +C.C., 

where A is the (real) amplitude, k is the wave number, and the temporal frequency ui is related to 
uq and k via the dispersion relation 

(3.3) u! = f'{uo)k — k^. 
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We now wish to describe the modulation of the linear wave packets. That is, we consider 
solutions of of the form (j3.2p . where the parameters A, k, and a; are allowed to depend on 

space and time. First, note that if A, k, and uj are allowed to vary significantly over one wavelength 
of the “carrier” wave (j3.2p . taken with A, k, and u constant, then there is no reason to expect 
that the evolving wave packet should be described as a modulated periodic wave. Hence, we only 
consider modulations, where A, k, and uj vary slowly on space and time. This naturally leads one 
to consider two separate space and time scales, namely the “fast” oscillations of the carrier wave 
and the slow evolution of A, k, and uj. 

Utilizing a multiple scales approach, we thus seek a modulated plane wave solution of (13.ip of 
the form 

(3.4) v{x, t, X, S) = A{X, ^ 

where (x, t) denote the “fast” variables and {X, S) = (ex, et) denote the “slow” variables. Here and 
elsewhere, e refers to a small, but non-zero, quantity. Substituting this into (j3.1h and collecting 
powers of e we find that 

(3.5) {-iuj + ikf'{uo) - ik^) + + (^o - “ ('“o “ Sk'^t))^ 



The 0(1) equation simply requires that the slowly evolving temporal and spatial frequencies satisfy 
the dispersion relation (j3.3l) . i.e. 

uj{k) = f{uo)k-k^, k = k(X,S), 

for all X,S G M. Since this forces uj to be real valued regardless of the slow evolution of the 
wave number k, we see that the modulated solutions (|3.4I) are bounded for all time, and hence the 
constant solution uq is considered to be modulationally stable. 

To describe the slow evolution more precisely, note that at 0(e) the expansion (13.51) yields two 
equations for the slow dynamics of k and A: 

/o dk ,. dk „ dA 

(3-6) ^ ^ 

where uj'(k) = /'(uq) — 3k‘^ depends on X and S through the evolution of k. Together with the 
dispersion relation, the system (|3.6I1 describes, to leading order, the evolution of slowly modulated 
wave packets in ()3.1I) . In particular, since uj(k) is real-valued for all A; G M, the linearization of this 
system about (k, A) = (/cq, Hq), corresponding a linear wave train of the form (|3.2n with k = k^ and 
A = Aq constant, is hyperbolic, and hence both k and A will remain bounded for all time. Again, 
this implies the modulational stability of the constant solution. 

Note that the Whitham modulation equations (j3.6p for slow modulations of a constant state 
prescribe the evolution of only two quantities, namely the wave number and the amplitude of the 
wave. Moreover, the evolution equation for k is independent of the evolution of the amplitude. 
These observations are due to the fact that the dispersion relation 13.31 directly relates the temporal 
and spatial frequencies. 

3.2 Nonlinear Dispersive Waves 

We now wish to understand the long-time dynamics of an arbitrary amplitude, slowly modulated 
periodic traveling waves of (j2.ip . The main technical difference is that the spatial and temporal 





frequencies are no longer coupled through the linaer dispersion relation, and hence we expect to 
find a modulation system describing the slow evolution of three of the characteristic parameters of 
the wave. 

Recall from Section [2] that periodic traveling waves of (12.ip are parametrized (up to transla¬ 
tions) by three parameters a, E, and c, and the period depends on the parameters. While this 
parametrization is natural from the existence standpoint, we will see below that it is is not well- 
suited to the formal asymptotic description of Whitham’s theory of modulations. (We ntilize this 
parametrization in the rigorous theory in the following sections, however.) Rather, the modulation 
equations are most naturally described in terms of the slow evolution of the period and partic¬ 
ular conserved quantities associated with the modulated solution. The Jacobian between these 
parametrizations will play a vital role in our theory. 

We seek an asymptotic description of modulated periodic traveling waves of m- To account 
for the separate fast and slow scales of space and time involved in such modulations, we utilize a 
multiple scales expansion. We introduce “slow” variables {X, S) = {ex, et) and note that, in the 
slow coordinates, (j2.1h becomes 

(3.7) Us = euxxx + f {u)x- 
Following Whitham [701172], we seek a solution of (|3.7p of the form 

(3.8) u{X, S) = ^\ x, 5 ) + + ©(e^), 

where the phase (j) is chosen to ensure that j = 0,1, are 1-periodic functions of the variable 
y = 4>{X,S)/£. In particular, we note that the local period of oscillation of the modulated wave 
train is £/dx(l>: where we assume the unknown phase a priori satishes that dx(j> 7 ^ 0. Notice 
that periodic traveling waves of (I2.1jl fit this asymptotic description by taking u = constant 
in the slow variables and choosing (j){X, S) = k{X — cS). That is, periodic traveling waves of (12.ip 
correspond to constant solutions in the slow coordinates. 

We now substitute (|3.8I1 into (|3.7I1 and collect like powers of e, resulting in a hierarchy of 
equations that must be satisfied. At the order of 0(e“^), we find that 

4>sdyU^^'> = {4>xdyfu^°'> + {(l)xdy)f{u^°'>) 

which is simply ()2.3[) for differentiated in the variable y, under the identification k = (f>x as 
the spatial frequency and c = u/k as the wave speed, where io = —(^>3 is the temporal frequency of 
the modulation, all of which are considered functions of the slow variables X and S. For hxed X 
and S, thus we may choose to be a periodic traveling wave of (12.ip . and hence of the form 

(3.9) A, 5) = u(y; a(X, S), E{X, S), c{X, S)) 

for some solution u of (|2.3p . where a, E, and c are independent of y. In particular, the compatibility 
condition cpsx = 4 >xs implies that 

(3.10) ks+ujx = 0, 

which serves as a nonlinear dispersion relation. It effectively replaces the dispersion relation for 
linear waves. The linear case satisfy (|3.10l) thanks to the dispersion relation and the first equation 
in (13.6h . In the literature, (|3.inp is sometimes referred to as the equation of conservation of waves. 
At the order of 0{1), continuing, 

(3.11) - dxf{u^^^) - 5x(</>^52)n(°), 
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supplemented with 27r-periodic boundary conditions, where 


= d\ — 2 : = ky, 

is recognized as the linear operator obtained from linearizing (I2.ip about the carrier wave 
The Fredholm alternative implies that (|3.1ip has a solution, provided that the right hand side is 
orthogonal to the kernel of the adjoint operator {dyC[u^^'^])^ = —C[u^^'^]dy in the sense of L^(0, 27r). 
Since is not constant, it can readily be shown by differentiating (12.31) that 

ker(£[ti^°^]9y) = span{l, 

Hence (13.lip has a solution, provided that the orthogonality conditions 

fMs + %G(«(»>)=0, 

' 1 Fs + 8xQ(«<»>) = 0 


hold, where 


M{X,S)=[ uo{y,X,S)dy and P{X,S)=[ ul{y,X,S)dy 
Jo Jo 

are functions of the slow variables X and S, and the fluxes G and Q denote the inner product 
of the right hand side of (13.lip with 1 and uq, respectively; in particular, they are ?/-independent 
functions of X and S that may be written out explicitly, if desired; see m, for instance. 

Together, three equations in ()3.10l) and ()3.12p form the Whitham averaged system for (|2.ip . 
describing the mean behavior of the slowly varying functions a, E, and c in p.9p implicitly in terms 
of the slow evolution of the functions k, M, and P. In particular, it formally follows that (|2.ip has 
a modulated wave train solution that to leading order can be described as 


u{x) = a(X, S),E{X, S),c{X, S)^ + 0(e), 

where a, E, and c evolve in the slow variables in such a way that k, M, and P evolve according to 
the (nonlinear) first order system (I3.inp - (l3.12p . 

To ascertain the modulational stability of periodic traveling wave uo{ao, Eq,co) of (12.ip . we seek 
to understand whether the functions a, E, and c in (j3.9p evolve away from [uq, Eo,cq) in the slow 
scale or not. The modulation system (I3.10p - (|3.12l) . however, describes the evolution of a, E and 
c implicitly through the evolution of k, M, and P. As such, if we want to measure deviation of 
(a, E, c) from (oq, Eq, cq) in terms of the evolution of {k, M, P), we must require that, at least in the 
vicinity of the carrier wave uq, nearby periodic traveling waves of (j2.ip can be re-parametrized by 
the parameters k, M, and P. By the implicit function theorem, such a (local) re-parametrization 
is possible, provided that the Jacobian matrix 


(3.13) 


d{k,M, P) 
d{a, E, c) 


is non-singular at {uq, Eq,co). 

Under the assumption that the matrix (I3.13p is non-singular, the modulational stability can 
be analyzed by linearizing the modulation equations (I3.10l) - (j3.12p about the constant solution 
(ao, Eq, Co). We seek solutions of (j3.10p - (j3.12p of the form 


(a(X,S),E(X,S),c(X,S)) 


(ao, Eo, Co) + e{d{X, S), E(X, S),c{X, S)), 


10 

































which, upon substituting into (I3.10p - (|3.12p and dropping all O(e^) term^, yields the linear system 

(d{kQ,Mo,PQ) , 9(a;o, Go, Qo) ^ ^ r, 

I 8(a,E,c) Bs(a,E.c)+ a„(a.E,c)) = 0, 

where the 0 subscript denotes evaluation at (oq, Eq,cq). This is a constant coefficient, linear system 
of PDE’s, and we may seek solutions of the form 


{d,E,c){X,S) = e 


XS-ikX 


(ao,Eo,co), A G C, kG 


for some (ao,Eo)Co) G M^, in which case the temporal and frequencies of the perturbation can be 
determined from the dispersion relation 


V„(X,k) := det ( p(^o.Mo.-Po) _.^d{u,„ao,Qo) 


d{a, E, c) 


d{a, E, c) 


= 0 , 


a homogeneous cubic polynomial in the variables ^ and k. In particular, k) = /c//i) 

and hence if {zj}^^^ are the roots of 

(3.15) Vw{z):=Vw{l,z) 

then the dispersion relation defines three spectral curves 

Xj{k) = izjk, j = 1, 2, 3. 

In particular, uo{ao, Eq,co) should be modulationally stable, provided that the reduced dispersion 
relation has three real roots, while uq should be modualtionally unstable if it admits at least one 
complex root with non-zero imaginary part. 

The above formal analysis suggests that a necessary condition for the modulational stability 
of a periodic wave train of p2.1l) is for the linearized Whitham system p3.14l) to be hyperbolic, 

i.e. the matrix has all real eigenvalues. In the next section, we rigorously derive a 

^(^0) ^0; ^o) 

modulational stability theory whose main result, Theorem 14.21 relates the stability of a periodic 
wave train of (12.ip to slow modulations to the eigenvalues of a 3 x 3 matrix. These eigenvalues end 
up agreeing precisely with the roots of providing a rigorous justification of the above formal 
calculations; see Theorem oi below. 


4 Rigorous Theory of Modulational Instability 

We discuss how to make rigorous the formal arguments sketched in the previous section in various 
settings. In the process, we provide a general methodology of detecting modulational instabilities 
for a large class of nonlinear dispersive equations. 

4.1 Analytic Setup 

Let no(-; a, E, c) denotes a T = T{a, E, c) periodic traveling wave solution of the generalized gKdV 
equation (|2.1I) . The idea of “stability” can formally be stated as requiring that if a solution to (|2.ip 
starts close to the solution uq then it stays close to the solution uq for all time. To understand 

^Note the 0{1) terms cancel out due to the definition of (ao, Eq, cq), so the leading order term is 0{e). 
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the dynamics of solutions that start near uq, one typically begins by linearizing the governing PDE 
about the solution uq, a process we now describe. 

First, notice that uq is a stationary T-periodic solution of the PDE 

(4.1) Ut = CUz + Uzzz + fiu)z, 

which is simply the gKdV equation (12.ip written in the traveling spatial frame z = x — ct. We now 
seek a solution u of with initial data u{0 ) = uq + vq for some vq G H^(M) with ||uo||j:f3(R) 
sufficiently smaljl. Writing v{t) = uq — u{t) for t > 0 for which u{t) is definecO, it follows that v 
satisfies the evolution equation 

(4.2) vt = Lv+Af[uo,v], 

where here L = duF{uo) : C L^(M) —>• L^(R) is the linear operator 

L = dz {dl + c + /'(uo)) , 

and where A/^Pwo, u] = dz (fiuo + v) — /(uq) — f'{uo)v) is a nonlinear operator satisfying the quadratic 
estimate 

for some constant C = C'(||'Uo||hi(r)) > 0 independent of v. We expect that as long as vq is 
sufficiently small in, say, L^(R) then the solution v{t) should be small at least for small, again 
in L^(R) say, at least for short times. Moreover, so long as v{t) remains sufficiently small it is 
reasonable to expect that the nonlinear evolution (j4.2p may be well approximated by the linear 
evolution equation 

(4.3) Vt = Lv 

posed on C L^(R). 

Due to linearity, we may of course relax the smallness assumption in the analysis of (j4.3p . 
The underlying wave uq is said to be linearly stable (to perturbations in L^(R)) provided that all 
solutions of (14.3p remain bounded in time, while it is linearly unstable (to L^(R) perturbations) if 
there exists some initial data uq G such that the associated solution of ()4.3I) grows in time. 

The linear stability of uq may be further reduced to a problem in spectral analysis as follows. 
As the linear evolution (14.31) is autonomous in time, we may take the Laplace transform in time, 
resulting in the problem 

(4.4) Lv = Xv, 

where here v{z) = v{z; A) denotes the Laplace transform of the function v{t, z) in (14.31) and A G C 
is the Laplace (temporal) frequency. Equation (|4.4p is a spectral problem posed on L^(R), and as 
such the stability of uq may be studied by analyzing the spectrum of L considered as a densely 
defined operator on L^(R). Indeed, if the L^(R)-spectrum of L intersects the open right half plane 
in C, then one expects the solution of (14.4p to grow in time and hence the underlying solution uq 

^Tiie requirement that the perturbation be integrable is not the only choice possible. Other natural candidates of 
classes of perturbations include periodic perturbations with fundamental period nT for some n = 1, 2, 3,.... While 
such periodic classes of perturbations are natural from a variations! viewpoint, they may impose artihcial constraints 
on the physical problem. As we will see below, the class of localized (i.e. integrable on the line) perturbations includes 
information about all quasi-periodic perturbations. 

^Here, we ignore issues of well-posedness of 1231 to such initial data. 
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should be unstable. With this in mind, we say that the underlying wave uq is spectrally unstable 
(to L^(]R) perturbations) if there exists a A with 3?(A) > 0 in the L^(M)- spectrum of L. Otherwise, 
we say uq is speetrally stable (to L^(M) perturbationsjl. Notice that since the spectral problem 
(j4.4ll is left invariant under the transformations 

V V and A i—>• A, 

where here the “bar” denotes complex conjugation, and 

z — 2 : and A —A 

the spectrum of L is symmetric with respect to reflections about both the real and imaginary axes. 
Consequently, rio is spectrally stable if and only if the L^(R)-spectrum of L is confined to the 
imaginary axis. 

As (|4.4p is a spectral problem for a linear operator with periodic coefficients, its spectrum 
is most conveniently characterized via Floquet-Bloch theory. Standard results from this theory 
[191 [62] dictate that non-trivial solutions of (|4.4p can not be integrable on R, i.e. they can not have 
finite norm in LP(R) for any 1 < p < 00 , and that they can at best be bounded functions on the 
line. Further, any bounded solution of (14.41) must be of the form 

v{x) = e''^^w{x) 

for some w € Lpgj.([0,T]) and ^ € [—tt/T, vr/T). In particular, a given A G C belongs to the L^(R)- 
spectrum of the operator L if and only if there is a non-trivial solution of the quasi-periodic spectral 
problem 

{ Lv = Au, 
v{x + T) = e*^'^u(x) 

for some ^ G [—'kIT^ttIT). Alternatively, A G C if and only if there exists a ^ G [—tt/T, vr/T] such 
that there is a non-trivial T-periodic solution of the equation 

L^w = \w, where := 

For each ^ G [—tt/T, vr/T), the spectrum of the operator L^, considered as a densely defined operator 
on Lpgj.([0, T]), is comprised merely of discrete eigenvalues with finite multiplicity and, furthermore, 

spec£2(]R)(L) = U specL2e^([o,T])(^s)- 

?e [-tt/T,tt/T) 

The parameter ^ is referred to as the Bloch frequency and the operators are called the Bloch 
operators associated to L. 

Rather than analyzing the essential spectrum of the operator L directly, the above discussion 
indicates that we may rather choose to study the T-periodic eigenvalues of the Bloch operators for 
each ^ G [—vr/T, vr/T). From a practical level, one does not expect to be able to explicitly compute 
the eigenvalues of for an arbitrary Bloch frequency except outside of a few very special case^. 
Thankfully, however,for the purposes of modulational stability analysis we need only consider the 
spectrum of the operators in a neighborhood of the origin in the spectral plane, and only for 

^Generally, spectral stability does not imply linear stability, as is familiar from basic ODE theory. Nevertheless, 
spectral instability often does imply linear (and nonlinear) instability. 

®For example, completely integrable systems m- 
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1^1 <C 1. Indeed, notice that if the eigenvalues of Lq are confined to the imaginary axis, then uq is 
spectrally stable to perturbations which are themselves T-periodic. When 0 < |^| <C 1, the period 
of the perturbation is neary that of the underlying wave uq, and hence the effect of the perturbation 


is seen only on very large space and time scales. Hence, the spectrum of the operators with 
0 < 1^1 <C 1 describes the stability of uq to long wavelength perturbations. 

Modulational perturbations clearly fall into the catagory of long-wavelength perturbations, but 
they form a special subclass where the effect of the perturbation is simply to modulate the defining 
characteristics of the wave, here given by the parameters a, E, c, and the translational mode zq. 
The key point is that variations in these defining parameters provide spectral information of Lq at 
A = 0 in the spectral plane. Formally, this can be seen by differentiating the profile equation (12.2p 
with respect to a, E, c, and z and noting that, 

(4.5) LQspsxi{dzUQ,daUQ,dEUQ} = 0, LqOcUq = -d^UQ 

where here the T-periodic boundary conditions have not yet been enforced; see Lemma l4. II below 
for a precise description of the generalized kernel of Lq. 

From the above remarks, it is natural that the spectral stability of the underlying wave uq to 
slow modulations corresponds to ensuring the spectrum of the Bloch operators L^ near A = 0 is 
confined to the imaginary axis for all 0 < |.^| <C 1. Our program is thus two-fold: first study the 
structure of the generalized kernel of the unmodulated operator Lq, and then use bifurcation theory 
to examine how the spectrum near the origin of the modulated operator Lg bifurcates from A = 0 
for 0 < 1^1 <C 1. This program is carried out in the next section. 

4.2 Modulational Instability in gKdV Type Equations 

We begin our program by demonstrating that, under a generic non-degeneracy condition, the 
unmodulated operator Lq has a two-dimensional T-periodic kernel with a Jordan chain of height 
one. Henceforth, we employ the notation 




and 


for determinants of 2 x 2 and 3x3 Jacobians, respectively. 

Lemma 4.1 ([ T 2 l ITSjL Suppose that u{x] oq, Eq, cq) is a T-periodic solution of the profile equa¬ 
tion (ESI) and that the Jacobian determinants Te, {T,M}a,E and {T, M, P}a,E,c are non-zero at 
{oq^Ecq). Then the functions 


f'o = {T,'a}a,E 
= {T,M}a,EUx 
(j)2 = {u,T,M}a,E,c 


V'O = 1 

V'l = Jq fi2(s)ds 

= {T,M}a,E + {T,M}a,EU 


are all T-periodic and satisfy the relations 


Lq4>q — LqcPi — 0 
Lq4>2 = —fii 



In particular, the functions form a basis for the generalized null-space of Lq, and the 

functions {V’i}j=o form a basis for the generalized null-space of L^. Moreover, ([ot]) ~ ^ 

ifiy^j- 
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The proof follows from a straightforward application of the Fredholm alternative; see [mils], 
for example. The main observation is that while Ux is T-periodic, the derivatives Ua, ue, and Uc 
generically are not T-periodic due to the dependence of the period on the parameters {a,E,c). 
Indeed, the change in these functions across a period is proportional to derivatives of the period: 
for example, 


/ ue{T) \ 

UEx{T) 

'^Exxi'^) 

\ ••• / 


( ue{^) \ 

UEx{^) 

V ••• / 


( ue{0) \ 
UExiO) 
UExx ( 0 ) 

V ••• 


with similar expressions for the change in Ua and Uc across a period. It is then trivial to check that 
the functions T-periodic and, thanks to (14.5p . satisfy the desired relations under the 

action of Lq. Further, since (j )2 clearly has mean-zero, the functions clearly T-periodic 

and, by again differentiating the profile equation, are seen to satisfy the desired relations under the 
action of Tg. Since u can be chosen, through appropriate translation, to be an even function, and 
since Lq is parity reversing, we can further verify that {ipj, ^2 qq -r]) = 0 for all j, fe = 0,1, 2 

which, by the Fredholm alternative, clearly implies the desired structure of the generalized kernels 
of To and lJ. 

Henceforth, we assume that the quantities Te, {T, M}a,E, and {T, M, P}a,E,c are non-zero at 
the underlying T-periodic traveling wave uq. Lemma 14.11 then implies that the elements of the 
T-periodic kernel of the unmodulated operator Lq are given by the elements of the tangent space 
to the (two-dimensional) manifold of solutions of fixed period and fixed wave speed at the point 
uq, while the generalized T-periodic kernel of Lq consists also of a vector in the tangent space 
to the (three-dimensional) manifold of solutions of fixed period at uq with no restrictions on the 
wave speed. It immediately follows that the origin is a T-periodic eigenvalue of Lq with algebraic 
multiplicity three and geometric multiplicity two. 

We now turn to study how this triple eigenvalues bifurcates from the (A,^) = (0,0) state. To 
this end, we note that a standard Baker-Campbell-Hausdorff expansion reveals that 


L^ — Lq + if,Li -|- L 2 + (if,)^ 

Lq = dxC is the unmodulated operator, Li := [Lq,x\ = C — 2df., ancH L 2 := [Li,x] = —3dx- Note 
that these operators are well-defined on Lpgj.([0, T]) even though the function x is not. Furthermore, 
a very useful observation is that 


(4.6) ~ ^ 

whenever i -|-j -|- fc = 0 mod(2), which can readily be verified by using Lemma l4.ll and noting that 
Lq and L 2 reverse parity while Li preserves parit}j§. Since L^ is a relatively compact perturbation 
of Lq depending analytically on the Bloch frequency it follows that the operator L^ will have 
three eigenvalues {Aj(^)}^^^, defined for |^| <C 1, bifurcating from A = 0 for 0 < |^| <C 1. To 
determine whether the underlying wave u is modulationally stable or not, it remains to track these 
eigenvalues as functions of ^ and, in particular, determine whether these eigenvalues are confined 
to the imaginary axis for all 0 < |^| <C 1 or not. 

^That is, Lo sends odd functions to even functions, and even functions to odd functions, all while preserving 
periodicity. 

®Here and throughout, [A, B] = AB — BA denotes the commutator of operators A and B. 

®That is, Lo and L 2 send even (odd) functions to odd (even) functions, while Li sends even (odd) functions to 
even (odd) functions. 
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Notice that by Lemma HTTl we may use the functions and {'ipj}‘j=o to construct explicit 

rank 3 eigenprojections 


no : -^^per([0> T]) gker(Lo), Hq : T]) gker(L|j) 


onto the total eigenspac^i^ for the operators Lq and Lg. We can thus represent the action of the 
operator Lq on generalized kernel of Lq as the 3x3 matrix operatoi0 


/ 0 0 0 

Mq := noLoIIo = I 0 0 {'ijji, Lq4>2) 
\ 0 0 0 


Notice that (V’l, Lo02) = ^{T, M}a^E{T, M, P}a,E,c / 0 by assumption, reflecting the Jordan struc¬ 
ture of the unmodulated operator Lq. Further, the above projections are non-degenerate in the 
sense the action of the identity operator 

/ {^o,4>o) 0 0 \ 

noHo = 0 0 

\ 0 0 0 ^ 2 , 4 > 2 ) ) 

is non-singular with determinant —{T,M}‘^^{T,M,P}a,E,c- 

Since A = 0 is an isolated eigenvalue of Lq, it follows by standard spectral perturbation theory 
that the above eigenprojections can be continued for |^| 1 into analytically varying rank 3 

projections 

3 3 ^ ^ 

: Ll,M T]) ^ 0 ker (L^ - A, (?)!), T]) ^ 0 ker 

i=i i=i 

with ranges coinciding with the total left and right eigenspaces associated with the eigenvalues 
{-^i(0}^=i of [5I]- fo particular, we may find analytically varying bases {uj(^)}^^g and 

of the total right and left eigenspaces for the eigenvalues in such a way that 

Uj(0) = cj)j and Uj(0) = 'ijjj. Note in particular that we need not require the {ujj^^g and {vj}‘j^Q 
to be eigen-bases for the total eigenspace, a degree of flexibility that is helpful during explicit 
calculations. Using these ^-dependent projections, for each fixed |^| <C 1 we can project the infinite 
dimensional spectral problem (14.4p onto the three dimensional total eigenspace associated with the 
{-^i ( 01^=1 in particular, the action of the operators L^ on this subspace can be represented 
by the 3x3 matrix ooeratoJ^ 

;= rt^L^rt^. 

The key observation is that for each fixed is a constant matrix and that, further, the eigenval¬ 
ues of the matrix coincide precisely with the eigenvalues for the modulated operators 

L^. In particular, for each |^| <C 1 the eigenvalues {Aj(.^)}^^;^ are the roots of the cubic characteristic 
equation 

det — An^rtg^ = 0. 

^°Here, the total eigenspaces for Lo and LJ are precisely gkerfLo) = ker(Lo) and gker = ker(Lo), respectively. 
^^Specifically, Mo := Lo(f)j) >^2 ([o,T])j and similarly for the projection of the identity; see below. In 

Section 15. 4. 21 we will utilize a slightly different rescaling of these matrix operators. 

^^As above, Mf := \{vjiP, LfVj{P) 1^2 froTid and similarly with the projection of the identity operator. 
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Note that since noHo is invertible, it follows by continuity that II^II^ is non-singular for all |^| <C 1 
. Thus, the three eigenvalues bifurcating from the triple eigenavlue A = 0 of Lq can be studied for 
for 1^1 <C 1 by analyzing the eigenvalues of the 3x3 matrix Mg, which is a seemingly much easier 
task. 

The eigenvalues of the matrix Mg may now be studied using standard techniques in matrix 
perturbation theory. First, the matrices Mg depend analytically on ^ for |^| <C 1, due to the ana- 
lyticity of the associated eigenprejections Ilg and Ilg and the analytic dependence of the operators 
Lg on and can be expanded for |^| <C 1 as 

Mg = Mo + (iOM(i) + (ie)^M(2) + OdeP). 

for some 3x3 matrices M^^) and M^^) whose coefficients are independent of It is well known 
that the eigenvalues of Mg in general bifurcate continuously in the perturbation parameter but 
may not vary in a manner. Rather, due to the Jordan block at ^ = 0, these eigenvalues are 
generically expected to admit a Puiseaux series expansion in fractional powers of In order to 
guarantee that the eigenvalues are at least in ^ near (^ = 0, it is sufficient to verify that the 1,2 
and 3,2 entries on the matrix M^^^ both vanish. By a straightforward calculation, it can be shown 
that 

= (^0, ^i</>i)L2^^([o,r]) and [M^%,2 = {ip2, Lih) > 

both of which vanish thanks to (I4.6p . Consequently, the eigenvalues are functions of 

^ and hence admit an expansion of the form 


(4.7) XjiC) = 1^1 < 1, 

for some function determined by the roots of the characteristic equation 


(4.8) 


det ( T^Mg — /rllgllg 


= 0 


for 1^1 <C 1. It follows that, to leading order, the modulational stability of the underlying wave may 
be inferred from determining the values of Hj{0) G C in (|4.7p . Note that if we further assume the 
numbers /Uj(0) are distinct (which generically holds in our case), then the eigenvalues {Aj(^)}|^]^ 
depend analytically on ^ for |^| <C 1. 

To determine the values in (14.7p . notice the Dunford integral calculus implies that the 

bases and have the same regularity on ^ as the eigenvalues and hence 

they may be expanded near ^ = 0 as 

(4.9) Vj{C) = 4>j + + o(6> VjiC) = ipj + + o(0 

for some functions and vp'^ independent of We wish to use these expansions to compute 
asymptotic expansions of the projections Ilg and Ilg valid near ^ = 0. Note, however, the non-trivial 
Jordan structure of Mq implies an inherent degeneracy in the limit ^ 0. Indeed, note that while 
the kernel of Mq is two-dimensional, there are three eigenvalues of defined near the origin for 

^ = 0. Moreover, the three eigenvectors of Mg tend to the same as ^ 0. An effecient way 

of unfolding this degeneracy is to rescale the matrices and Ilgllg as 

Mg := iE(0-'MgE(0, /g := 

^®This is a consequence of the fact that all three of the eigenfunctions for associated to the eigenvalues {Aj(^)}|^i 
tend to (/>! as 5 —>■ 0, a fact seen both analytically and numerically. 
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where 


0 0 

S(0 := I 0 1 0 
0 0 zC 

By straightforward computations, it follows that for |^| <C 1 we have the expansions 
^ {iI)q,Li(I)q) mi,2 {ll)Q,Li(t)2) + {v^q\Lq(I)2\ \ 


Mf = 


+ L(iv\ 


( 1 ) 


{'lpl,Lo4>2) 




+ 0(1) 


where 


mj,2 := {i^j,L2ct)i +Liz;f^^ + (vf\Li(t)i i = 1,3 


and that, similarly. 


I, = 


^ i'lpoAo) + (v^o\(pi^ 0 ^ 

0 0 
^ 0 (i;2,v^^'''^ + (v^^\4>i'^ {i^2,4>2) J 


+ 0 ( 1 ). 


Therefore, it follows that the numbers /ij(0) are determined as the roots of the (cubic) effective 
dispersion relation 

V{ff) := det ^Mq — /z/q^ = 0 

This establishes the following theorem, which is a restatement of Theorem 2 in [45]. 

Theorem 4.2 (|l5])' Let zz(-; oq, T^o, cq) be a T = T{ao, Eq,co) periodic traveling wave solution of 
(|2.1[) and suppose that Te, {T,M}a,E, and {T, M, P}a,E,c are non-zero at oq, Eq, and cq. Then 
for |,f| <C 1, the triple eigenvalue A = 0 of Lq bifurcates into 3 eigenvalues Xjif) that are in ^ 
near ^ = 0 and satisfy 

XjiO = iTjC + o{C), j = 1,2,3 

where the constants fij € C are the eigenvalues of the effective dispersion matrix 


(4.10) 


D(a, E, c) 



{a,E,c) 


at {a,E,c) = {ao, Eq,co). In particular, a necessary condition for u to he modulationally stable is 
that the dispersion matrix be weakly hyperbolic, i.e. has three real eigenvalues. 

To utilize Theorem 14.21 in practice, one must determine the first-order correctors z;g^\ and 
in ()4.9p . We note, however, that in the small amplitude limit these correctors do not contribute 
to leading order and hence are not necessary to compute the roots of V\ see Section 15.4.21 below. 
Outside of this asymptotic regime, however, these correctors play an important role. First, notice 
that the functions {z;j(^)}y^Q and {z}y(^)}y^Q satisfy the relations 

Bg (Lgu,(0) = hvffO, ng (lJu,(o) = L\vffO, 

respectively, for all |^| <C 1. In particular, expanding in ^ we find that 

+ (-^1 + ni+o) 4>j € Range (11(0)) 
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and 


Llvf^ + [l\ + HiLy V'j G Range (n(O)j 

In particular, since = 0 we may choose so that = —Li4>i. An explicit formula for 

can be found by noting that, by definition, 

(4.11) LiUx = [Lo,x]ux = Lq {xUx) 

and that TexUx + Tue is T-periodic. Similarly, since ker(LQ) = span{V'j}j=o ,2 we may choose 
and such that = ~L\'tpj for j = 0,2. Unlike explicit forms of j = 0,2 

are not needed; their dehning relations are sufficient for calculation purposes. With these choices 
one can, through tedious but straightforward calculations, compute the effective dispersion relation 
explicitly, yielding 

T^il^) = C + I ({T, P}e,c + 2{M, P}a,E) - 1{T, M, P}a,E,?j 

for some non-zero constant C = C{a,E,c). The reality of the roots of T> follows from the sign of 
its discriminant, given here as 

Ami ■■= I ({T, P}e,c + 2{M, PU,Ef - M, P}Ie,c- 

In particular, the polynomial P has three distinct real roots provided that Ami > 0, while it has 
one real and two (non-real) complex conjugate roots if Ami < 0. It follows that if Ami < 0, 
then uq is modulationally unstable with one branch of spectrum bifurcating from the origin along 
the imaginary axis two branches bifurcating in complex directions. It follows that a necessary 
condition for modulational stability is Ami P 0. Further, when Ami > 0, a simple symmetry 
argument implies that the three spectral branches bifurcating from the origin are confined (locally) 
to the imaginary axis, corresponding to a triple covering, yielding modulational stability. 

Corollary 4.3 ([12]). Under the hypotheses of Theorem If.Sj a necessary condition for the mod¬ 
ulational of the underlying wave u Ami > 0. Furthermore, a sufficient condition for u to be 
modulationally stable is that Ami > 0. 

4.3 Connection to Whitham Modnlation Theory 

In Section [3] we used formal asymptotics to determine a dispersion relation Vw in (I3.15P with 
the property the linearized Whitham system (I3.14[) has eigenvalues of the form A = ikp whenever 
/i G C satisfies Vw{z) = 0. In Section above, we demonstrated that he eigenvalues for the Bloch 
operators Lg expand as A = ifif, + o(^) where ^ G C is one of the three (distinct by assumption) 
eigenvalues of the effective dispersion matrix 'D{a,E,c) defined in (j4.10p . These two approaches 
are in fact equivalent for equations of gKdV type. 

Theorem 4.4 ([3HI35]). Under the assumptions of Theorem ff.Sl the eigenvalues of the effective 
dispersion matrix (I4.10p agree with the roots of the reduced Whitham dispersion relation D]v{z) = 
0. In particular, a necessary condition for the modulational stability of the underlying periodic 
wave train is that the linearized Whitham system (I3.14h be weakly hyperbolic, or, equivalently, that 
equation Vw{z) = 0 has three real roots. Further, a sufficient condition for modulational stability 
is for (I3.14P to be strictly hyperbolic, i.e. the equation Vw{z) = 0 has three distinct real roots. 
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Theorem 14.41 rigorously justifies formal predictions of Whitham’s theory of modulations in the 
context of gKdV equations of the form (12.11) . It follows that a sufficient condition for the modula- 
tional instability of the unerlying carrier wave is the ellipticity of the associated linearized Whitham 
system. Its proof follows by a long and tedious calculation demonstrating that there exists a con¬ 
stant C / 0 such that 

det (D(a, E, c) — /rl) = CVwifJ-) 

for all /i € C: see [441 145] for details. We merely pause to note that, in a sense, Theorem 14.41 
implies that the processes of averaging and linearizing commute: the linearized Whitham system 
being obtained from and averaging process (Fredholm alternative) and then linearizing, while the 
rigorous theory obtained from linearizing the PDF and then applying the Fredholm alternative 
(averaging). This observation suggests the above program can be readily extended to equations 
outside the context of ( 12 .ip 

One class of equations for which the above program has been particularly successful is in the 
context of dissipative PDF of the form 

Ut +f{u)x + g{u) = {B{u)ua;)^, X E M, u E M” 

consisting of a continuum of models, with / = 0 corresponding to systems of reaction diffusion 
equations, and (7 = 0 to systems of conservation laws. Not only has Whitham’s theory of modula¬ 
tions been rigorously validated in this context at the level of spectral stability, it has also provided 
key insights into the nonlinear dynamics of modulated waves in such models, as well as the stability 
and dynamics of wave trains to nonlocalized perturbations, consisting of a localized part plus an 
asymptotic phase shift at the spatial infinities; see [m n sg sa 146] . These additional insights 
come from continuing the asymptotic WKB expansion in Section [3| to next order in e, obtaining a 
diffusive correction to the Whitham modulation system. The details of this program are beyond 
the scope of the current manuscript, and the reader is referred to the above references for details. 


4.4 Evaluation of A 


MI 


In Section 4.2 above, we defined a modulational instability index Ami, the sign of which determines 
the modulational stability of a given periodic traveling wave solution of the gKdV equation (|2.ip . 
To compute Ami in practice, we must determine derivatives of the quantities T, M, and P with 
respect to the traveling wave parameters a, E, and c, a task that in general seems formidable. 
Nevertheless, when / is a polynomial nonlinearity these derivatives may be computed explicitly in 
terms of the underlying wave. Indeed, note that if / is a polynomial then the quantities T, M, and 
P can be expressed in terms of the moments 


Cfc := f u''{x)dx = (f 
Jo J-y 


u 


yjE -v{u-a,c) 


du, 


where here E — V{u\a,c) is a polynomial of one degree higher than /; indeed, T = 70 , M = 71 , 
and P = 72 - Further, the derivatives of T, M, and P may be expressed in terms of the moments 


Ik := 


7 {E-V{u-,a,c)f^ 


du; 


for example, Te = lo and Pc = ip- The amazing fact is that if / is a polynomial of degree n, 
then the 2 n — 1 moments expressed as linear combinations of the n — 1 moments 

{CfcjfcZo’ ^ known as the Picard-Fuchs relation. As such, the modulational instability index 
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Am/ can be explicitly computed in terms of moments of the underlying wave u; see m for more 
details. 

To illustrate how this process works, suppose / is a polynomial of degree n — 1 and set P{u) = 
E — V{u; a, c). Then we can write P{u) = oq + aiu + 02 //^ + ... Onu” for some constants 
Defining the moments Cfc and 1^ as above, it clearly follows that 

dCk ^ ^ ^ j 
dttj dok 2 

for each k,j = 0,1, 2,... , n. Upon observing that 


Cfc — 


u^P{u) 

P{uf/^ 


n 


du = 2_^ o.iI. 
j=o 


m = 0,1, 2,... , n — 1 


and tha10 


2fcCfc-i = 2k 


,k-l 




du = 


u^P'{u) 

P{uf/^ 


du = 

j=0 


j+k — l: 


k = 0, 


, n — 1 


we find the following linear system of 2n — 1 equations in the unknowns {Ik}‘^Q^' 


(4.12) 
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Tldn 

0 
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hn-S 





V 

0 


0 

Ol 

202 


nan / 


\ l2n-2 

/ 


2(n - l)Cn-2 

/ 


which we refer to as the Picard-Fuchs system. The matrix that arises in the above linear system 
may be recognized as the Sylvester matrix of the polynomials P{u) and P'{u), A standard result 
in commutative algebra states that the Sylvester matrix of two polynomials is non-singular if and 
only if the two polynomials have no common roots. Thus, the above matrix is invertible provided 
that the polynomial P has only simple roots, in which case the moments {Ik}^^'^ can be explicitly 
expressed as a linear combination of the {Ck}^Zo with coefficients explicitly determined by the 
coefficients of P. 

For a given polynomial P, the above linear system can be solved easily with the aid of a computer 
algebra system. In the next section, we discussion applications of the above procedure in the case 
where P{u) = E — V{u; a, c) and V is the effective potential energy associated to (12.111 for various 
nonlinearities /. 


5 Applications 

In Section 4 above, we defined a modulational instability index Am/, the sign of which determines 
the modulational stability of a given periodic traveling wave solution of the gKdV equation (12.111 . 
We now discuss a few examples where one can use the Picard-Fuchs system (j4.12l) to explicitly 
evaluate Ami in terms of the underlying wave. We note that these examples are not exhaustive 
but are chosen for illustrative purposes. 

^^Here, we take the convention that the left hand side is zero when k = 0. 
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5.1 The KdV Equation 

Consider first the KdV equation 

(5.1) Ut — Uxxx {u )x- 


In this case, using notation from Section O the effective potential energy is 

c 1 

V(u; a, c) = —au — -u^ + -v? 

^ O 

and the periodic traveling waves u{x‘,a,E,c) of (|5.ip exist provided that the discriminant 

disc {E — V{-;a, c)) = — (l6o^ + 3a^c^ — 36Eac — 6Ec^ — 36E^) 

is Dositivj^. Further, the only periodic solutions along the surface disc (K — V(-; a, c)) = 0 are 
constants. In this case, setting 


Cfc := 


du 


\/2{E — V(u; a, c) 


h = 


du 


{2{E-Viu-,a,c))f^' 


we find that (T, M, P) = (Co, Ci, C 2 ) and that 


^ E,a,cCk — 


Ifc, lfc+1, 2^1 k+2 


In particular, the Picard-Fuchs system (|4.12p becomes 


/E a § 1 0 \ 


/ Io\ 


/ T \ 

0 K a § I 


h 


M 

a c —1 0 0 


h 

= 

0 

0 a c —1 0 


h 


2T 

\ 0 0 a c -1 / 


\ h y 


V 4M / 


which, after some elementary matrix algebra, implies the following formulae for 


Te 


(4a + c^)M + (6K + ac)T ^ [M jP-q.c) 

12disc(K-V(-;a,c)) ’ ^ 12disc (K - V(-; a, c)) 


{r,M,p} 


a^E^c 


T^{E -V{M/T-,a,c)) 
2disc (E — V (•; a, c)) 


and 


2Ami 


+ a2,iT‘^M + ai^2TM‘^ + 
2i037disc(K-V(-;a,c))^ 


see m- Here the coefficients Oij are given by 


“ 3,0 = 36K + ISaEc — 8a^, 02,1 = 18Ec^ — 6a^c + 36aE, 

02,1 = —18cE + 24a^ + 3ac^, 00,3 = + 6ac + 12£'. 

^®Positivity of the discriminant ensures the polynomial E — V(u;a,c) has three real roots in u. 
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The monotonicity of the period, Te > 0, holds by a result of Schaaf[63], while {T, M, P}a,E,c > 0 
clearly holds. Further, {T,M}a,E is positive due to Jensen’s inequality and the fact that 

/ V'{u; a, c) du ^ 
yjE — V{u-,a,c) 

Thus, the hypothesis of Theorem 14.21 are satisfied, and hence Corollary 14.31 implies that the sign 
of Ami determines the modulational stability of any underlying wave u. Since the numerator of 
Ami is seen numerically not to vanish, it follows that Ami > 0. Taken together, this yields the 
following. 

Theorem 5.1 f|13jl. All periodic traveling wave solutions of the KdV equation dSI]) are modula- 
tionally stable. 


5.2 The Modified KdV Equation 


Next, we consider the modified KdV (mKdV) equation 
(5.2) Ut — Uxxx T (p ')x^ 


which corresponds to the focusing (+) and defocusing (-) modified KdV (mKdV) equations. Like 
the KdV equation, (j5.2p is completely integrable and hence its sepctrum can, in principle, be 
explicitly computed with the aid of completely integrability and algebro-geometric techniques. Of 
course, as mentioned above, our methods are independent of integrability. Further, as we will see 
our methods yield a modulational instability result for (15.211 that seems to have been previously 
missed. 

As discussed in Section [2] above, the existence of periodic solutions for (|5.2p differs between the 
focusing and defocusing cases. In each case, however, existence may be inferred from properties of 
the associated effective potential energy 


14 (u;a,c) := =f-m^ 



— au 


Indeed, in the focusing case, periodic solutions exist provided that the quartic polynomial E — 
V+(-;a,c) has at least two real roots and if its discriminant is non-zero. In the defocusing case, 
however, periodic solutions only exist when the discriminant oi E — a, c) is positive, indicative 
of the polynomial having four distinct real roots. See Figure [Din Section [2l 
In this case, the Picard-Fuchs system reads 
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Tl 
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\ h / 
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6P 

/ 


which can be readily solved to provide explicit formulas for Te, {T,M}a,E-, and {T,M,P}a,E,c in 
terms of T and P anc0 the parameters {a,E,c)-, see m- Further, we find 

» ^mKdv{T,P,a,E,c)‘^ 

^MI = - o 

disc {E — V±(-; a, c)) 

^®The fact that M does not enter the formulas for the gradients of T, M, P with respect to a, E, c is easily explained 
by complex analysis; see m for details. 
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where VmKdV'^^ a homogeneous polynomial of degree 6 in T and P with coefficients depending on a, 
E, and c. The sign of this clearly agrees with the sign of the discriminant of the quartic polynomial 
E — V±(-; a, c), which in turn is positive if it has four distinct real roots, and negative if it has only 
two distinct real roots and two distinct complex conjugate roots. This leads us to the surprisingly 
simple characterization of modulational instability for the mKdV equation. 

Theorem 5.2 ([13] i. The traveling wave solutions to the niKdV equation (1521) are modulationally 
unstable for a given set of parameters (a, E, c) if the polynomial 

c 1 

V{u-, a,E,c) = E + au + -u^ ± 

has only two distinct real roots and two distinct complex conjugate roots, and it is modulationally 
stable if it has four distinct real roots. 

As mentioned above, the structure of the roots of V{--, a, E, c) depends greatly on whether the 
equation is focusing or defocusing. In the case of the focusing mKdV, Theorem 15.21 implies that if 
the parameter values give rise to only one periodic solution, then this solution is modulationally 
unstable. Conversely, if the parameters give rise to two periodic solutions, then both of the solutions 
are modulationally stable. In terms of the explicit solutions described in Section [2. II above, we have 
the following corollary. 

Corollary 5.3. For the focusing mKdV equation, all the dnoidal wave solutions, constructed in 
Section{2A\ above, are modulationally stable, while all the cnoidal waves are modulationally unstable. 

In contrast, notice that for periodic waves to exist in the defocusing mKdV V must have four 
distinct zeroes. This simple observation yields the following. 

Corollary 5.4. All periodic traveling wave solutions of the defocusing mKdV equation are modu¬ 
lationally stable. 

In particular, Corollary 15.41 implies the snoidal wave solutions described in Section 12.11 above 
are modulationally stable. 

5.3 The Schamel Equation 

Lastly, we study the Schamel equation (after taking u >->■ —u in ()5.1h i 
(5.3) Ut + Uxxx + 2 ~ 

where u must be taken to be positive. Periodic traveling waves of (15.31) are determined implicitly 
via 

/ du 

® —= X — ct. 

J2{E + art + 

Using the obvious change of variables u = v‘^ shows that the traveling wave solution to this equation 
may in turn be found from 

l^ = E + av^ + -v^-v\ 

2 2 

Thus, the physically admissible parameter regime corresponds to the set of all (a, E, c) G for 
which there exists a bounded interval in (0, oo) such that the quintic polynomial E + av'^ + — v^ 

is non-negative, corresponding to a non-negative periodic solution. 
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In this case, the Picard-Fuchs system is the set of 9 equations 
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where 

f 2v^ dv 

Ck ■= f —? A: = 0,1,2,... 

Y 2{E + av‘^ + \cv^ — v^) 

and Ik are defined similarly; see m, for instance. In particular, notice that T = ^i, M = ^3 and 
P = Cs- While the expressions are too long to list here, they are easily handled with the aid of a 
computer algebra system, in which case one finds that Ami > 0 in the entire domain of existence. 
In conclusion, all periodic traveling waves of (15.31) are modulationally stable', see m for details. 


5.4 Extensions to Equations with Nonlocal Dispersion 

The analysis in previous sections readily extends to other equations allowing for nonlocal dispersion. 
We will illustrate this by discussing equations of KdV type 

(5.4) ut = Mux +f{u)x, 
where A4 is defined via the Fourier series 

u{x) = => (A4n)(x) = ^ m(n)h(n)e*"'®, 

characterizing dispersion in the limit limit. We assume throughout that the symbol m is even and 
real-valued. 

In the case of m{k) = —k"^, notably, ()5.4p recovers the generalized KdV equation. Notice, 
however, that (j5.4p is nonlocal unless the dispersion symbol m is a polynomial of ik. Examples 
include the KdV equation with fractional dispersion 

(5.5) ut = + f{u)^, -l<a>2, 

where A = ^ —df. is defined via its Fourier multiplier m{n) = |u|. In the case of a = 2, it recovers the 
generalized KdV equation, and in the case of a = 1 and f{u) = , it corresponds to the Benjamin- 

Ono equation, which was proposed in [3] and m as a model of the unidirectional propagation of 
internal waves of small amplitudes in deep water. Recalling that the Fourier multiplier of the 
Hilbert transform H is —isgia{k), the Benjamin-Ono equation may alternatively be written 

(5.6) ut = Hd^u + uux- 

Like the KdV and modified KdV equations, the Benjamin-Ono equation is completely integrable, 
although, of course, our methods do not rely on this fact. In the case of a = —1/2 and f{u) = 
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moreover, (15.41) was argued in [38] to have relevance to water waves in two dimensions and in 
the infinite depths. Incidentally fractional powers of the Laplacian occurs in many applications, 
such as dislocation dynamics in crystals (see [T8|, for instance) and financial mathematics (see 
m, for instance). Notice that the dispersion symbol of (15. 5 p is homogeneous in the sense that 
m{Xk) = a^m{k) for all A; G M and a > 0, which makes the analysis for waves of all amplitudes 
tractable by the methods described in the previous sections. 

However, many applications require non-homogeneous symbols for their dispersive operators. 
The intermediate long-wave (ILW) equation 


(5.7) 


Ut + Ux + {l/H)Ux - NhUx + {u^)^ = 0. 


is such an example, where here 77 > 0 is a parameter and the Fourier multiplier Mh is defined via 
its symbol m{k] H) = kcoih.{kH). The ILW equation is indeed of the form (15.4p . where the symbol 
is given by m(^) = 1 + 1/H — Mh-, and it arises as a model to describe the interface between density 
stratified fluids, both with finite depth [50] . 


Another equation with non-homgeneous dispersion corresponds to the choices m{k) := y 
and f{u) = n^, in which case (15.4p becomes the Whitham equation. The Whitham equation was 
originally put forward by Whitham in [721 p. 477] to model unidirectional propagation of surface 
water waves with small amplitudes, but not necessarily long wavelengths, in a channel. Notice that 


iti«i 

so that we may regard the KdV equation (15.ip as to approximate up to “second” order the dis¬ 
persion relation for the Whitham equation, and hence for the full water wave problem, in the 
long-wavelength regime (i.e. small frequencies h). Unlike the KdV, however, the Whitham equa¬ 
tion may be used to model short wave uni-directional phenomena. As a matter of fact, Whitham 
advocated that this equation could be used to explain wave breaking and peaking, both of which 
are unobserved in the KdV approximation. In light of the famous Benjamin-Feir instability, it 
is natural to expect that small periodic wave trains of the Whitham equation are modulationally 
unstable provided their wavelength is sufficiently short. 

Unlike the gKdV equation (12.111 . the existence of periodic traveling waves for nonlocal equations 
of the form (j5.4l) is not trivial. Indeed, a traveling wave of (15.4p is a solution u{x, t) = u{x — ct) 
where u{-) satisfies, after one integration, the profile equation 


(5.8) 


— Aiu -|- cu -|- f{u) = a 


for some constant a G M. For a broad range of dispersion symbols and nonlinearities, a plethora of 
periodic traveling waves of (|5.4p may be attained from variational arguments, e.g. the mountain 
pass theorem applied to a suitable functional whose critical points satisfy (15.811 . In each example 
considered below, there will exist a smooth four-parameter (including translations) family of peri¬ 
odic traveling waves, in agreement with the existence theory in the local gKdV case considered in 
previous sections. 

Once existence and an appropriate parametrization of the periodic solutions of ()5.8I) has been 
established, modulational stability can be studied through studying the L^(M)-spectrum of the 
linearized operator 

L := dx (-M + c - /'(n)) 
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near the origin. As in the case of local gKdV equations, since the function u is T-periodic, the 
L^(M)-spectruni of L is readily described by Floquet-Bloch theory: 

spec^2(]u)(-t') = U 

^e[-7r/T,7r/T) 


the Bloch operators := act as usual on Lpgj.([0, T]). The modulational stability theory 

developed in Section 0] can thus be applied directly in this case provided the Jordan structure of 
the generalized kernel of L is analogous to that described in Lemma l4.ll 


5.4.1 The Benjamin-Ono Equation 

Let us consider the Benjamin-Ono equation (after the change of variables u —u in (IS.Gp ) 
(5.9) ut - Aux + {u^)x = 0. 


In [3], Benjamin exploited the Poisson summation formula and found that for each A; > 0 there 
exists a two-parameter (up to translations) family of 27r/A;-periodic traveling waves of (15.9p of the 
form 


(5.10) 


u{z; a, k, c) = , k -^(y/c^ -4a -h c), 

c2-4a ^ ^ 

:2-4a-F - 


where a and c are arbitrary constants constrained by the conditions 


(5.11) 


c < 0 and k‘^ < — 4a. 


It follows that (j5.9p admits a three parameter family (up to translations) of periodic traveling wave 
solutions, and that these solutions can be parametrized by the period and the parameters a and c. 
Here the period, equivalently k, is a parameter that is independent of a and c so that, in particular, 
derivatives of the solution with respect to a and c are automatically T-periodic in x. Since dSjD 
obeys Galilean invariance under 

(5.12) u{z; a — cX s‘^,k,c — 2s) = u{z; a, k,c) X 


for all A G M, upon an appropriate choice of A, one may assume that a = 0. 

To analyze the modulational instability of a given solution u of the above form, we turn our 
attention to the Lpgj,([0,T]) spectrum near the origin of the Bloch operators 

Lg := e-^^^dx (A + c - 2u) 

for 1^1 <C 1. Although the symbol of A is not smooth near the origin, the symbol of dxA is at 
the origin and we can expand the Bloch operators as 

= Lq {if)Li L2 + o(I^P) 


where 

(5.13) To := dx {A-2u- c) =: dxC 
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and 


Li = [Lq,x] = 2K — 2u —c, L2 = [Li,x\=K ^dx- 

Following the analysis of Section [H the first step in modulational stability analysis is to understand 
the Jordan structure of the kernel of the unmodulated operator Lq. This is achieved in the following 
lemma. 

Lemma 5.5 ([H] ). Suppose that u = u(-; oq, feo, cq) is a 2tt /k-periodic traveling wave solution 
of (15.9|] and that {M,P}a,c ■= det is non-zero at {ao,ko,co). Then the generalized 

Lpej.([0,T]) kernel of the linear operator Lq defined in (|5.13l) possesses the following Jordan block 
structure: 

(i) dim(ker(Lo)) = 2 
(a) dim(ker(Lo))/dim(ker(Lo)) = 1 

(Hi) dim(ker(Lg'’'^))/dim(ker(LQ)) = 0 for all integers n >2. 

In particular, the functions 

Vo = Ua Wo = McU - Pc 

Vi =Ux Wi= {MaUc - McUa) dx 

V2 = Uc W2 = Pa - MaU 

form a basis and dual basis for the generalized kernel of L, respectively: 


LqVo = Lvi = 0 , LqV2 = V2 

L\wo = P^W2 = 0 , L\wi = W2. 

Moreover, {wj,Vk) = —{M,P}a,c^jk and {wj,LiVk) = 0 whenever j + k-\-l is even. 

The proof follows the same lines as the proof of Lemma 14.11 see HU. The key step to extend 
the proof to the current nonlocal case is to verify that the kernel of Lq is only two dimensional. 
This, would follow essentially immediately provided one could verify that the operator £ in (|5.13p 
is non-degenerate, i.e. that 

(5.14) ker (£) = spanjua,}. 

While this is obvious in the local gKdV case, by classical Sturm-Liouville theory, it is far from obvi¬ 
ous in this more general nonlocal context. Nevertheless, by relating via complex analysis techniques 
the nonlocal profile equation (15.8|) to a fully nonlinear ODE, Amick and Tolland [1] demonstrated 
that the operator C associated to the Benjamin-Ono equation is indeed non-degenerate in both the 
solitary and periodic settings. Equipped with this non-degeneracy result, the proof of Lemma [52] 
follows from the same arguments as Lemma l4. II 

To verify the hypothesis of Lemma 15.51 owing to the explicit solution formulae (I5.inh , we note 
that a straightforward calculation (see m for details) reveals that 

f27r/k ^ - 

M(a, k, c) := / uiz; a, k, c) dz = 2tt — —(v — 4a J- c) 

Jo k 
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and 


r27vjk ^ 

P{a,k,c) :=- / a, k, c) dz = —ctt + jr(\/— 4a + c)^. 

2 Jq Ak 

so that, in particular, 


{M,P},,,= 


27r^ 


> 0 


A:-\/c?Tr^ 

for all {a,k,c) satisfying the constraints (I5.11|) . Lemma 15.51 thus provides an explicit basis for the 
generalized kernels of the unmodulated operator Lq and its adjoint about any periodic traveling 
wave solution of (|5.9p . 

With Lemma [5.51 in hand, the methodology of Section [4.21 mav be applied without modification. 
In particular, the triple eigenvalues at A = 0 of Lq bifurcates into 3 eigenvalues Aj(^) that are 
in ^ near ^ = 0 and satisfy 

XjiO = itijC + o{C), j = 1,2,3 

where the constants G C are the eigenvalues of, taking into account the parametrization of 
solutions in ()5.10p . the effective dispersion matrix D(o, A:,c) in ()4.10l) . Using (I5.10p this matrix 
can be calcnlated explicitly in terms of the traveling wave parameters (a, k, c) provided an explicit 
formula for the corrector in (14.911 can be found: due to the parametrization of solutions of (|5.9p 
by their period, the formula from our gKdV analysis mnst be modified. Since (j4.1ip still holds in 
the present case, may be identified as an appropriate multiple of the 211 /k periodic function 
xUx + Tut where T := see 1141 L emma 9]. With this choice, the matrix D(a, /c, c) can be shown, 
up to a similarity transformatiorj^d to be given (for a = 0) by 

—ttT (vtT)^ — (vr/c)^ 0 

D(a = 0, fc, c) = I 1 ttT 0 1 , T = 

27r^ 0 vrT 


2it 

T' 


A quick calculation shows the eigenvalues of D(0, fc, c) are irT and ztvrT \J2 — {cT')~‘^. Since periodic 
traveling waves of (j5.9p exist only when (fi > k = l/T, it follows that the eigenvalues of D(0, A:,c) 
are real and distinct. By the Galilean invariance (j5.12p and Theorem 14.21 this yields the following 
result. 


Theorem 5.6 ([E]). All periodic traveling wave solutions of the Benjamin-Ono equation dSj]) are 
modulationally stable. 

In fact, the above calcnlations on the Benjamin-Ono equation (15.9|] is merely a special case of the 
recent analysis of Bronski and Hur where they considered the modulational instability of arbitrary 
amplitude periodic traveling wave solutions to the fractional KdV equation (|5.5p with 1/3 < a < 2. 
For such a, one can show through the calculus of variations that for each T > 0 there exists a 
two-parameter (np to translations) family of periodic traveling wave solutions parametrized by the 
constant of integration a and the wave speed c and that, furthermore, these waves depend on a, T, 
and c in a manner. A key difficulty in this more general setting is to verify the non-degeneracy 
property (15.141) for all 1/3 < a < 2. Although the methods of Amik & Tolland [T], described above, 
do not apply outside the Benjamin-Ono case a = 1, it is possible to verify through a periodic 
adaptation of the recent nonlocal Sturm-Liouville theory of Frank & Lenzmann [29] for fractional 
Schrodinger operators to verify that the kernel of Lq is indeed two dimensional for all 1/3 < a < 2; 
see [HSQ] for details. Equipped with this non-degeneracy result, the effective dispersion matrix 
can be numerically computed for all such a; see [14] for details. 

^^Indeed, this formula for D was derived in [T4| using a direct spectral perturbation expansion, a method equivalent 
to that described in Section sn above. 
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5.4.2 The Whitham Equation for Water Waves 

As described above, in the absence of explicit solution formulas it is in general very difficult to calcu¬ 
late the effective dispersion matrix (I4.10p explicitly, and one must typically rely on well-conditioned 
numerical techniques. One case where this difficulty is avoided is in the small amplitude limit. In¬ 
deed, as we will see below, for a wide class of nonlocal evolution equations of the form (j5.4p such 
small amplitude solutions can be expanded (in amplitude) explicitly to arbitrarily high algebraic 
order. Equipped with such an explicit expansion of solutions, one can compute an aymptotic ex¬ 
pansion of the now amplitude dependent effective dispersion matrix (jd.lOp and hence ascertain the 
modulational stability for waves of sufficiently small amplitudes. 

To illustrate this procedure, we now consider the modulational stability of periodic traveling 
wave solutions of the Whitham equation 

(5.15) ut +Mux +{u^)^ = Q 

where here the symbol of the operator A4 is given by m{n) := described above, 

equation (I5.15P combines the full uni-directional dispersion relation for surface water waves with a 
canonical shallow water nonlinearity. Details of the forthcoming analysis can be found in [39]. 

Periodic traveling wave solutions of (|5.15p take the form u{x,t) = u{x — ct), where c > 0 and u 
satisfies 

(5.16) M.U — cu + u^ = {1 — c)‘^b 

for some constant 6 E M; the normalizing factor of (1 — c)^ is added for convenience. To describe the 
small amplitude periodic wave trains of (I5.15P we seek solutions of (I5.16P of the form u{x) = w{z) 
where z = kx, where here w is required to be 27r-periodic and /c > 0 is interpreted as a wave 
number. The 27r-periodic function w is then required to be a solution of 

(5.17) M-kW — cw + = {\ — c)^6 

where here is defined via 

Mke'"'^ = m{kn)e^"'^, n € Z. 

Observe that the symbol here is not homogeneous and hence the wavenumber may not be factored 
out of the symbol. 

Noting that is a continuous operator for all A; > 0 and s > 0, a simple 

argument verifies that solutions of (j5.17p with ||rt;||Loo < | belong to the Sobolev space and 
hence are smooth. As a result, smooth 27r-periodic solutions of (I5.17P with small amplitudes may 
be sought as roots of the function 

F{w, k, c, b) = M-kW — cw + up' — {1 — c)‘^b] 

note that F : x M_|_ x M_|_ x M ^ by a Sobolev inequality. Clearly F{0; k, c, 0) =0 for every 

k,c ^ M+ and, furthermore, the derivative of this map at such a trivial solution is the operator 

9^E(0,/c,c,0) = Mk -c 

which has a trivial kernel in provided that c ^ m{kn) for any n € Z. By the implicit function 
theorem, near any such trivial solution for which du,F has a trivial kernel, there exists only other 
trivial solutions. Thus, to construct non-trivial solutions of (I5.17P we start at a zero solution for 
which the kernel of d^F is non-trivial. 
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To this end, notice that if c = m{k) then 


ker {du,F{0, k, c, 0)) = span . 


Using a straightforward Lyapnnov-Schmidt reduction, it follows that for each k > 0 there exists a 
two-parameter (up to translations) family of non-trivial roots of F, corresponding to smooth (even) 
27r-periodic solutions of (j5.17p . Furthermore, their small amplitude asymptotics may be found by 
seeking solutions of (|5.17p of the form 


w{k, A, b){z) = WQ{k, b) + Acos( 2 ;) + A^W 2 {z, b) + A^w^{z, b) + 0{\A\‘^) 
c{k, A, b) = m{k) + A^C 2 {k, b) + 0{\A\^). 


where here ^ is a real parameter (a sort of generalized amplitude), the wj are taken to be even, 
and Qb is defined for |6| <C 1 via F{wo{k,b), k,m{k),b) = 0. Plugging the above expansions into 
(j5.17p we find 


w{k, A, b){z) = Wo{k, b) + Acos( 2 :) -|- -A 


1 


+ 


cos(2z) 


m{k) — 1) m{k) — m{2k) 


(5.18) 


+ 0 ( 71(712 + 62 )) 


c{k, A, b) = co{k, b) + A 


1 


+ 


m{k) — 1 2{m{k) — m{2k)) 

+ 0{A{A^ + b^)) 


as 1^1,16| —)• 0, where 


co(A:, b) = m{k) + 26(1 — m{k)) — 662(1 — m{k)) + 0(6^) 


and 

WQ{k, h) = 6(1 — m{k)) — b^{l — m{k)) + 0(6^); 

although higher order expressions are obtainable, they are not necessary in our analysis. Together 
with translation invariance, this yields a four parameter family of small amplitude periodic traveling 
waves of the Whitham equation (IS.lSp : see [39l[23l[23] for more details. Throughout the remainder, 
we take 6 = 0 for simplicity. 

With the above parametrization in hand, we now turn to study the +2,^ := Lperii^^ ^vr]) spectrum 
near the origin of the Bloch operators 

L^{k, A) := e-^^^L{k, 7l)e*«^ L{k, A) := d, {-Mk + c{k, 7l, 0) - 2w{k, 7l, 0)) 


for |(^,^)| <C 1. Unlike the case of the Benjamin-Ono equation considered in Section [5.4.11 above, 
the symbol m(-) is smooth near the origin and hence we may easily expand the Bloch operators as 


L^{k, 71) = Lo{k, 71) + {iOLiik, 71) + {i0^L2{k, T) + 0(|eP) 


where 


Lo(^) A) = dz {—Aik + m{k)) — 2adz (cos(z)-) + 0{A? + 62) 
Li(fe, ^4) = [—dzAik^ z] + 'm{k) — 2acos(z) + 0{A‘^ + b"^) 
L2{k, A) = [[—dzAik, 2 ^], ^] + 0{A‘^ + 62 ) 
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for |A| <C 1. As in our previous analysis, the key step in the modulational stability analysis is to 
understand the Jordan structure of the kernel of the unmodulated operator Lo(A:, A). To study the 
spectrum of near the origin, we first note that the operator estimate 

\\L^{k,A)-L^{kM\a^^ =0{A) 

27r 27r 

as A —>• 0, uniformly in ^ G [—1/2,1/2), which may be verified by bruit force, implies that the 
spectrum of L^(k,A) is close to that of the constant coefficient operator L^{k,0), uniformly in the 
Bloch frequency Consequently, we may consider the spectrum of L^{k, A) as a small perturbation 
of that of L^{k,A). In particular, the spectral projections and from Section [J] are now 
amplitude dependent, described to leading order by the projections associated to the trivial solution 
A = 0. 

We concentrate hrst on the right eigenprojection. Being constant coefficient, the spectrum of 
Lg(A:,0) for each k > 0 and ^ G [—1/2, 1/2) can be computed via Fourier analysis as 

spec 0)) = {iujn,^{k) : n G Z} C iM 

where ujn,^ik) := (u+^) {m{k) — m{kn + k^)) with corresponding eigenfunctions In particular, 
when ^ = 0 we have w_i^o = ^ 0,0 = <^ 1,0 = 0 so that A = 0 is an eigenvalue of L^i^k, 0) of algebraic 
multiplicity three, with eigenfunctions cos(z), sin(z), and 1. Further, it follows from above that 
the eigenspaces associated to the three non-zero eigenvalues w±i,g,Ci;o,g are independent of As 
such, the associated right spectral projection 11^(A:, 0) is independent of Finally, we note that the 
adjoint operators L^{k, 0), being constant coefficient, share the same eigenvalues and eigenfunctions, 
and hence the left spectral projection n^(A;,0) is independent of ^ as well. In particular, we finci^ 

Mo(A:, 0) := (n^L^n^) (A;, 0) = I —— w_i^g) -|-w_i^g) 0 1 . 

\ 0 0 / 

For |A| <C 1, it follows that the spectral projections can be expanded as 

U^{k, A) = no(A:, A) + 0(A2 + f), YL^{k, A) = no(A:, A) + ©(A^ + f) 

so that, to first order, the spectral projections for the left and right eigenspaces of L^{k,A) agree 
with those of the unmodulated operator LQ{k,A). In particular, the variations in these projections 
in ^ do not effect the leading order asymptotics and, furthermore, the bases used to construct 
these projections agree to first order. As such, to construct the projections n^(A:, A) and n^(/c. A) 
to leading order, we need only find a basis for the generalized kernel of Lo(A:,A) The origin is 
an eigenvalue of Lo(A:,A, 6) of algebraic multiplicity three and geometric multiplicity two for all 

^®Following [39] . the projections here are slightly rescaled from that of Section ID Here we use := 
/ {tpi, (jjj) ];^2 ([o,Ti)j and similarly for the projection of the identity. 
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1 ^ 1 , | 6 | <C 1 and, furthermore, the functions 

1 


4>i{z) := 


2(1 — m{k)) 


[{dbc)daw - {dac)dbw] {k,A,b){z) 


= cos{z) H- ^ ~ Qcos{z)b + 0{A‘^ + 6 ^), 


<(' 2 (^) := - -dzw{k,A,b){z) 
a 


= sin(z) + 


m{k) — m{2k) 
i,b){z) 
sin(22; 


A + 0(A^ + b^) 


m{k) — m{2k) 

(t)3{z) :=1 

form a basis of the generalized eigenspace for LQ{k,A,b) for all |yl|, | 6 | <C 1 ; specifically 

= Lo4’2 = 0 , and Lo </>3 = — 2002 . 

Using this basis, we find 


/ 0 

0 

0 \ / 

—km'{k) 

0 

0 \ 

M^(fe,A) := (n^L^n^ (fe,A)= 0 

0 

2 + id 

0 

—km'{k) 

0 

^ ^ \0 

0 

0 / \ 

0 

0 

m{k) — 1 / 


i^A( 


m{k) — 1 


V 2{m{k) — m{2k)) 


+ ( km'{k) + ^k‘^m''{k) 


0 0 2 
0 0 0 
1 0 0 

0-10 

1 0 0 1 +C)(|Cp + |^|^' 

0 0 0 


and, similarly, that 




A 


m(k)-m(2k) y ^ ^ 

Next, we turn to study the roots of the characteristic polynomial 
P(A, k, A) = det 


0 0 1 

0 0 0 I +0{A^) 


A) - A {k, A 

= C3(C, k, A)X^ + ic2{(-, k, A)X‘^ + ci(^; k, A)X + ico{C; k, A) 


for each fixed A: > 0 and (A,^, yl)| <C 1. As before, the roots of P{-,^]k, A) correspond to the 
eigenvalues of the Bloch operators L^{k,A) bifurcating from the origin for |(^, A)| <C 1. We first 
note that the coefficient functions cj depend smoothly on ^ and A for |(^, A)| <C 1 and that, 
moreover, the Cj are real-valued since the spectrum of L^{k,A) is symmetric about the imaginary 
axis. Further, the spectral symmetry 

speCi 2 ^(L 5 (A;, A)) = speCi 2 ^(L_^(/c, A)) 

implies that the functions C 3 and ci are even functions of ^ while C 2 and cq are odd, and that the 
Cj’s are all even in A. Since A = 0 is clearly a root of P{X, 0; k, A) with multiplicity three and ^ = 0 
is a root of P(0,^; A:, A) with multiplicity three for all |A| <C 1, it follows that 

Cj{(;k,A) = dj{C;k,A)f~^, j = 0,1,2,3 
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Figure 3: Plot of r(A:) for Whitham equation. 


for some real functions dj that depend smoothly on A and ^ and are even in A. Thus, the roots of 
P{-,C) may be written as A = —i^X where AT is a root of 


k, A) = det Al) + X {k, A)^ 

= k, A)X^ + d2{i; k, A)X^ + di{^-, k, A)X - do{^; k, A); 


compare to (I4.8p . It follows that the underlying periodic traveling wave is modulationally unstable if 
the characteristic polynomial P(— k, A) admits a pair of complex roots, i.e. if its discriminant 

^^,k,A (^^8d^d2didQ + d^d'i + T 4d3d^ — 27d3dQ) (i^, k, ^4) < 0 

for |(^, A)| <C 1, and that it is modulationally stable if the polynomial admits three real roots, i.e. 

if ^^,k,A > 0. 

Since the dj’s are even in A, we may expand 

^^,k,A = ^i,k,o + A{k)A? + 0{A?{A? + 


as |(^, ^)| <C 1, expressing that, to leading order, the modulational stability is governed by that of 
the limiting constant state at ^ = 0. Using the above explicit formulae, we find that = 0 

and 

(wo,.c - 

=- ^6 - > ^ 

for all 1^1 <C 1. Therefore, the sign of the discriminant A^^k,A for |(^,A)| <C 1 is determined by 
the sign of the function A(A:). Indeed, from above we see that if A(A:) > 0 then > 0 fro all 

|(^, A)| <C 1, implying modulational stability, while if A{k) < 0 then A^^^.A < 0 for |A| <C 1 fixed 
and 1^1 sufficiently small. The modulational stability of the underlying small amplitude waves is 
thus determined completely by the wavenumber k. 

With the aid of computer algebra software, we find 


where 


m 


{k{m{k) - iy{k{m{k) - 1)" ^ 

m{k) — m{2k) 


T{k) := 2{m{k) — m{2k)) + {k{m{k) — 1)'. 
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Recalling that m{k) := y , we see that the function k i-;- k{m{k) — 1) is strictly decreasing 

and concave down on the interval (0,oo), so that the sign of A(A;) is determined by the sign of 
r(fe). Further, a numerical evaluation of F shows that there exists a unique k* ~ 1.146 such that 
F(A:) > 0 for A: G (0, k*) and r(fe) < 0 for A; > A;*; see Figure [3l 

Theorem 5.7 ([39]). A 2t:/ k-periodic traveling wave solution of ()5.15p is sufficiently small am¬ 
plitude is modulationally unstable if k > k* k, 1.146, and modulationally stable if 0 < k < k*. 

Theorem 15.71 qualatativelv captures the famous Benjamin-Feir instability in the context of the 
Whitham equation: see |39| and |66] for more discussion on this connection. 


5.4.3 Fractional KdV and ILW Equations 

Finally, we note that the small amplitude analysis in Section [3] is largely independent of the specific 
symbol for the nonlocal operator A4. Indeed, the only features of the symbol m(-) that was used in 
the derivation of the modulational stability index A(A;) were (i) that m(0) = 1 and (ii) that k > 0 
is such that m{k) m{nk) for all n = 2,3,4,...; note that while (ii) provides a restriction on k, 
the condition m(0) = 1 can be relaxed via a simple scaling argumenl|l£j. 

Examples where this analysis can be applied without any modification include that fractional 
KdV (fKdV) equation, corresponding to (15.511 with f{u) = and the intermediate long-wave 
equation (|5.7I) . For the fKdV, we find for eaclo « > ^ that 


AfKdv(^; a) 


2A:4“a(l + a)^(2"+i - 3 - a) 
2 “ - 1 


which is negative (positive) for all A: > 0 provided that a < 1 (a > 1), indicating modulational 
instability (stability) of the associated small amplitude periodic wave trains. Recall that a = 1 
corresponds to the Benjamin-Ono equation, which was analyzed in SectionjATT] The corresponding 
result for more general power-law nonlinearities may be found in [43j . 

Similarly, for the ILW equation (15.71) we find for each H > 0 that 


^iLw{k] H) 




1) cosh(i7A:) -|- cosh(3iLA:) — 8Hk sinffiHk))'^ 
32H^ sinh(iLA;)i2 


^iLw{Hk) 


where Tjlw{z) := 1 — 2z^ — cosh(2z) -|- 2zsinh(2z); see [39|. A numerical evaluation indicates that 
^iLw{z) > 0 for all z > 0, and hence that, for each H > 0, the small amplitude periodic traveling 
wave solutions of (15^1) are modulationally stabl® 


6 Concluding Remarks 

We have made a detailed survey of many recent techniques rigorously studying the modulational 
instability of periodic traveling waves of KdV type equations. They are very general and apply 
to many other equations including systems involving dissipation. They are independent of the 
integrability of the governing equations. The theoretical approach in Section [4] may be used to 
rigorously validate formal predictions from Whitham’s theory of modulations in a variety of settings; 
see |9|, |58|, for instance, for more concrete examples, other than those listed in Section (4.31 We 

Although, the formula for A(fc) must be modified accordingly. 

^'^The condition q > 1/2 is an artefact of the corresponding existence theory. See [43] for details. 

^^This rigorously justifies a formal ’’amplitude equation” calculation in ED. 
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applied our theory to a number of concrete examples in Section [5l To illustrate the robustness 
of our theory, we discussed extensions to nonlocal equations of KdV type and obtained detailed 
results for a number of concrete examples. 

We have been concerned with the linear (spectral) stability of periodic traveling waves of KdV 
type equations to slow modulations. A fundamental open problem is to understand the nonlinear 
dynamics associated with modulationally stable waves and, even more interestingly, induced by a 
modulational instability. They have been studied in the context of dissipative systems (see [38] . 
for instance, and references therein), but it is not clear how such techniques apply to dispersive 
equations such as the gKdV equations. This is an exciting direction for future work. 
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